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At long distances interactions between neutral ground state atoms can be described by the Van der 
Waals potential V(r) = — X^^e C n /r n . In the ultra-cold regime atom-atom scattering is dominated 
by s- waves phase shifts given by an effective range expansion pcot So(p) — — 1/qo + rop /2 + . . . in 
terms of the scattering length cto and the effective range ro. We show that while for these potentials 
the scattering length cannot be predicted, the effective range is given by the universal low energy 
theorem ro = A + B/ao + C/a^ where A,B and C depend on the dispersion coefficients C n and 
the reduced di-atom mass. We confront this formula to about a hundred determinations of ro and 
ao and show why the result is dominated by the leading dispersion coefficient Ce ■ Universality and 
scaling extends much beyond naive dimensional analysis estimates. 
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Van der Waals (VdW) forces appear ubiquitously in 
many contexts of atomic, molecular, nuclear and particle 
physics. They account for long range dipolc fluctuations 
between charge neutral atomic and molecular systems [lj 
with implications on the production of Bose-Einstein con- 
densates of ultra-cold atoms and molecules ■ The in- 
termediate range nucleon-nucleon interaction due to two 
pion exchange also exhibits this VdW behaviour based 
on chiral symmetry Q providing a justification for the 
liquid drop model of nuclei Q . The short distance gluon 
exchange interaction between (colour neutral) hadrons 
also display this kind of interaction 0, [f| . Van der Waals 
forces, however, diverge when naively extrapolated to 
short distance scales @, Q. The study of such prob- 
lems in a variety of situations will certainly shed light on 
the usefulness of renormalization ideas within the specific 
context of quantum mechanics (see e.g. Ref. [§]). 

Fundamental work for neutral atoms was initiated in 
Refs. 0, 11, lH (see also 13 1), within a quantum-defect 
theoretical viewpoint. In this letter we systematically 
show that these simplified approaches work and analyze 
why they succeed. VdW forces are extremely simple in 
this case and are described by the potential 
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where C„ are the VdW coefficients which are computed 
ab initio from intensive electronic orbital atomic struc- 
ture calculations (see e.g. Ref. 14] for a compilation). 
Usually, only the terms with n = 6, 8, 10 are retained al- 
though the series is expected to diverge asymptotically, 
C n ~ n\ [15|| . The impressive calculation in Hydrogen up 
to C32 [l6| exhibits the behaviour C n ~ (l/2)™n! at rela- 
tively low n-values. The potential ([1]) holds for distances 
much larger than the ionization length l\ = h/ y/2m e I 
(I is the ionization potential) which usually is a few a.u. 
In the Born-Oppcnhcimcr approximation the quantum 



mechanical problem consists of solving the Schrodinger 
equation for the two atoms apart a distance r, 



< + U(r)u k 



1(1 + 1) 



Uk = k 2 Uk ■ 
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where U(r) = 2/iV(r)/fi 2 is the reduced potential, fi — 
m\mij(rn\ + m-i) the reduced di-atom mass, k = p/h = 
2tt/X the wavenumber, and Uk(r) the reduced wave func- 
tion. For our purposes, it is convenient to write the re- 
duced potential in VdW units, 



U(r) = - 



Rt 



92^ 
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where R e = (2^,C 6 /H 2 ) 4 is the VdW length scale and 



{2fxC 6+2n /h 2 )R 6 



for n > 1. In Table U we 



9n 

display numerical values for several di-atomic systems 
which are extremely small in the VdW units g\ ~ 10~ 2 
and gi ~ 10 -4 with the exception of H and He. Thus, 
we may anticipate a dominance of the Cq term in the 
calculations at low energies. This will be systematically 
quantified below. 

Using the superposition principle for positive energy 
scattering s- waves we decompose the general solution as 



u k {r) = u k , c (r) + k cot S (k)u k , s {r) , 



(4) 



with Uk, c ( r ) — > cos(fcr) and Uk, s (r) — > sin(fcr)//c for 
r — > 00 and So(k) the scattering phase shift for the I = 
angular momentum state. The potential given by Eq. (TT]) 
is both long range and singular at short distances. At 
short distances, the De Broglie wavelength is slowly vary- 
ing and hence a WKB approximation holds yielding 
for r — * 



Uk(r) 



Rn 



sm 



(5) 
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TABLE I: Van der Waals length Re = (2fiC 6 /h 2 )i, and the 
coefficients gi and 32 defined by the reduced dimensionless po- 
tential 2^V{r)Rl/h 2 = -x~ 6 [1 + gix" 2 
x — r/Re. We use results from Refs. [l7l 



with 



Atoms 


Re,(a.u.) 


.91 (10- 2 ) 


S2(10- 4 ) 


Li-Li 


64.9214 


1.42458 


2.97874 


Na-Na 


89.8620 


0.92320 


1.11369 


K-K 


128.9846 


0.64780 


0.49784 


Rb-Rb 


164.1528 


0.45647 


0.23370 


Cs-Cs 


201.8432 


0.36544 


0.13983 


Fr-Fr 


215.0006 


0.27362 


0.09526 


Li-Na 


73.2251 


1.25605 


2.17183 


Li-K 


84.2285 


1.18374 


1.79689 


Li-Rb 


88.0587 


1.18572 


1.70555 


Li-Cs 


92.8950 


1.21364 


1.68241 


Na-K 


106.5708 


0.80600 


0.80155 


Na-Rb 


115.3377 


0.74528 


0.65923 


Na-Cs 


123.2277 


0.73874 


0.61148 


K-Rb 


142.8292 


0.56543 


0.37106 


K-Cs 


154.2909 


0.53903 


0.32152 


Rb-Cs 


180.8480 


0.41520 


0.18654 


Be-Be 


43.3013 


2.54953 


6.70303 


Mg-Mg 


72.3589 


1.26414 


1.60399 


Ca-Ca 


111.4907 


0.84584 


0.65057 


Sr-Sr 


148.9023 


0.55117 


0.27632 


Ba-Ba 


189.4340 


0.41692 


0.15260 


Cr-Cr 


91.2731 


1.22821 




H-H 


10.4532 


17.51760 


423.45426 


He-He 


10.1610 


9.35937 


117.94642 



where R n — (2fiC n /h 2 ) 1 ^ n ~ 2 ^ corresponds to the highest 
VdW scale considered in Eq. ([1]). The phase (fk is in 
principle arbitrary. On the other hand, at low energies 
one has the effective range expansion |23j 



k cot So(k) 



1 



-r k 2 +v 2 k 4 \og{k 2 



(6) 



where ao is the scattering length, and ro is the effective 
range which can be calculated from 
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2 / dr 
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Here, the zero energy solution becomes, from Eq. (|4"|), 

uo(r) = u 0<c (r) - u 0<s (r)/a , (8) 
where Uo,c(?") — > 1 and UQ iS (r) — > r for r — + 00, yielding 



A B C 

a q.q 



(9) 



with A,B and C given by 



A = 2 dr(l - < c ) : 
Jo 



D 



dr(r - u , c uo,s) , 



C = 21 dr(r' - < s ) . 
'0 



(10) 

(11) 
(12) 



Then, combining the zero and finite energy wave func- 
tions we get for any r c > 

00 r°° 

u 'k( r ) u a{r) - u' a {r)u k {r) = k 2 / u k (r)u {r)dr ,(13) 

where r c plays the role of a short distance cut-off which 
is innocuous provided li <C r c <C R n - Using Eq. ((U, 
Eq. ^ and Eq. (jHJ) we then get for r c — > 



1 



1 f 

— sm((p k -(f ) = k 2 dr 

Rn Jo 

x Uk, c (r) + kcot5 (k) Uk, a (r) .(14) 



uo, c (r) u , s (r) 

a 



Orthogonality between Uk and uo requires (fk = <po in 
which case, expanding the integrand we get the structure 



kcot5o(k) 



a Q A{k)+B{k) 
a C(k)+V(k) 



(15) 



where A, B, C and T> are suitable functions of k. The 
interesting feature of Eqs. © and (fT5"|) is that the de- 
pendence on the scattering length ao and the potential 
is explicitly disentangled. This is a universal form of a 
low energy theorem, which applies to any potential reg- 
ular or singular at the origin which falls off faster than 
1/r 5 at large distances. We can visualize Eq. © as a 
long distance (VdW) correlation between r and ao. If 
the reduced potential depends on a single scale i?, i.e. 
U(r) = —F(r/R)/R 2 , one gets universal scaling relations 



r 4 = A + BK + cX 



R 



(16) 



where A,B and C are purely geometric numbers which 
depend solely on the functional form of the potential. 
For the pure VdW case, V = —Ce/r 6 , the effective 



range has been computed analytically [lOj, and in har- 



mony with the general structure Eq. (|16p reads 

v f? R^ 

-f = 1.39473- 1.33333— + 0.637318^ 

R e a ag 



(17) 



In Fig. Q] we confront the prediction for the effective range 
to the result of many potential calculations in VdW units. 
As can be vividly seen the agreement is rather impres- 
sive taking into account the simplicity of Eq. (fT7|) . Al- 
ternatively and discarding the exceptional outliers for ao 
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FIG. 1: The effective range ro vs the inverse scattering length 1/cto in units of the VdW radius 7?6 = for different 

ranges. Points are potential calculations [3, J2I, [1| (Li-Li,Na-Na), [H (Cs-Cs), [H (Na-Rb), [51 (Be-Be), [lj (Cs-Rb), UK 
(Cr-Cr), [13 (Fr-Fr), [H 0, HI Ell (H-H), (He-He). The line corresponds to Eg. ifTTll [i(llTll|. 
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FIG. 2: The effective range coefficients A, B and C (see 
Eq. Q) as a function of the dimensionless coupling gi repre- 
senting the 1/r 8 correction to the VdW potential 1/r 6 . 



and ro, we perform a \ 2 fit to the form Eq. (fT6| for the 
remaining 82 points and get A = 1.31, B = —1.57 and 
C = 0.66 in good agreement with Eq. (|17p . To our knowl- 
edge, Fig. [T] represents a completely unforeseen universal 
correlation supported by phenomenology. 

We analyze the robustness of this agreement by show- 
ing in Fig. [2] the effect on the effective range coefficients, 
A(gi), B(gi) and C(gi) due to adding a 1/r 8 term. From 
Table[T] we see that mostly g\ ~ 10~ 2 , in which case tiny 
changes are expected from Fig. [2] Actually, the small- 
ness of the deviations suggests using perturbation theory. 
If we expand the full solutions of the 1/r 6 potential at 
small k, Uk(r) — u (r) + /c 2 u 2 (r) + . . . , the change in 
the effective range due to inclusion of a AU(r) ~ 1/r 8 



potential keeping ao fixed reads 0, [37J 

/>oo 

Ar = 4 / AU(r)u (r)u 2 (r) dr 



(18) 



where r c ~ <?ii?6 "C i?6- The leading contribution is 
determined by the short distance behaviour of Uo(r), see 
Eq. ([5]), and we find i*2(r) ~ (r/R^) 4 uo(r) yielding Aro ~ 
g±Re log gi + 0(g±). The logarithmic enhancement in g\ 
can indeed be observed in Fig. [2] by the tiny curvature. 

The Cg-dominance is compelling and we show in Fig. [3] 
the universal functions A(k), B(k), C(k) and T>(k) which 
in conjunction with ao allow to determine the phase-shift 
from Eq. (fT5|) . They scale with Rq and are uniquely de- 
termined by the power law — 1 /r 6 once and forever. We 
have found that these functions show little dependence 
on gi and gi at momenta as large as fci?6 ~ 10, a rather 
unexpected result, hinting that the VdW universality and 
scaling extends much beyond the naive dimensional anal- 
ysis estimate fci?6 ~ 1 or the effective range expansion of 
Eq. (|6|) truncated with the first two terms. We note in 
passing that although such a truncation suggests a higher 
degree of universality, the VdW nature of the interaction 
prevents using (ao,ro) as fully independent variables in 
view of Eq. (|17|) and Fig. Q] It is interesting to notice 
that within a Nuclear Physics context characterized by 
short ranged Yukawa potentials stemming from meson 
exchange [38| the general Eq. © has been exploited as a 
means to check SU(4) Wigner symmetry in the Nucleon- 
Nucleon interaction for the 1 5'o and 3 S\ channels with a 
pattern similar to Fig. [T] 

We conclude by underlining that, when suitably dis- 
played, the analytical approach to Van der Waals forces 
pursued in previous works 1(J 11 1 acquires a quite uni- 
versal character with indisputable phenomenological suc- 
cess; the leading 1/r 6 contribution suffices to accurately 
describe low energy atom-atom scattering with just two 
parameters in a wide energy range. We naturally ex- 
pect new universality and scaling patterns to emerge 
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FIG. 3: The universal functions in units of denned by 
2fjV(r) = —Rg/r 6 which allow to determine the phase shift 
if the scattering length is also known, see Eq. (|15p . 

from systems characterized by power law forces but less 
understood such as molecular interactions in the ultra- 
cold regime. From a broader perspective we stress that 
the lack of dependence of potential model calculations 
on short distance details, unveiled from our system- 
atic comparison, is a feature traditionally built-in by 
the quantum-defect theory. Its natural counterpart of 
a smooth and controllable short distance cut-off depen- 
dence complies to the requirement of renormalizability 
within a pure quantum-mechanical framework. 
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